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BEHAVIOR OF SEPARATION CONSTANTS 
FOR FINITE GYROMAGNETIC PLASMAS 

Abstract 

! 

The type of the functions appearing in the expression for 
electromagnetic fields in finite magnetoplasmas depends only 
on the properties of the medium manifested in the separation 
constants. The co-)3 plane of the Brillouin diagram is 
partitioned into a definite number of zones characterized by 
the type of field solutions. The general quantitative behavior 
of the Separation constant for finite magnetoplasma in a longi- 
tudinal magnetic field on" the Brillouin diagram is described. 



BEHAVIOR OF GYROTROPIC PLASMA SEPARATION CONSTANTS 
AS FUNCTIONS IN THE co-j3 PLANE 



Of the numerous papers concerned with the propagation of waves 

through waveguides filled with gyrotropic plasma columns the majority 

2 3 4 5 

is confined either to the quasi-static approximation ' ' ' or to a dis- 

7 

cussion of the problem in general terms with tabulation of auxiliary functions , 

0 

such as, e.g. , the biradial functions and an outline of the general 

method of approach for obtaining the solution^ without a 

systematic and comprehensive presentation and classification of the 

numerous types of to “ propagation diagrams and field solutions such 

as are known, e.g. , for empty waveguides or waveguides filled with the 
11 12 

dielectrics and for slow wave structures. Only a few representative 
solutions have been published thus far^^'^"^'^^ and they did not include 
a complete systematic survey of all possible cases although they were 

I 

certainly oriented in this direction. 

Despite a high number of combinations that can arise with the 

several parameters that can be varied, viz. , the plasma density or the 

corresponding plasma frequency co , the intensity of the collimating 

P 

magnetic field or the corresponding cyclotron frequency , and the 

diameter of the waveguide or the corresponding empty waveguide- cutoff 
frequency, it is nevertheless possible to establish a system of classifi- 
cation of modes of propagation. 

It is believed that one way to make it possible readily to obtain the 
actual solutions for given conditions in waveguides is to prepare 



2 



catalogues where quantitites of interest such as, e.g. , the arguments of 



Bessel functions, biradial functions, and the like are tabulated and 
presented graphically on charts and diagrams as functions in the to- ^ 
plane for various combinations of values of the plasma frequency, 
cyclotron frequency, and the cutoff frequency of the empty waveguide. 

In addition to these, actual solutions for combinations of the above 
values in reasonable and suitably chosen steps would also be catalogued. 

The objective of this paper is to present a study of the separation 
constants, if indeed these quantities may be called by this name, which 
in circular cylindrical geometry appear as arguments of Bessel functions 

t 

in the expressions for the electromagnetic fields^ ' ^ viz. , 



E = 
z 









( 1 ) 



H = 



c^y(Ui)In(Ui y(u^)In(u3 ^ ) 

2 2 
^1 “ ^2 



( 2 ) 



Evidently, the separation constants determine the kind of the Bessel 
functions in Eqs. 1 and 2. If either u^ or u_ or both are imaginary, 

X Cm 

the corresponding Bessel function become modified, I instead of J . 

n n 

If they are complex conjugates Beseel functions of complex arguments 
appear in the above expressions. It is therefore believed that it is more 
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useful to know the behavior of the separation constants Uj^ and U 2 in 
the 0) - jS plane than merely the behavior of the separate components 
of the equivalent dielectric tensor c which is sometimes presented. 

The constants u^^ and u^ are solutions of a biquadratic equation 
which is obtained solely from the properties of the medium. In addition, 
of course, they must like, c^^ and c^, also be compatible with the 
boundary conditions. The medium or electronic equation was first 

g 

derived by Hahn and we shall give it in its original form: 






(3) 



The following notation is being used here: 




(4 -a) 





(4-b) 



k 

cy 




(4-c) 



(ji 



cy 





(4-d) 



k = 



(4-e) 



jW 

c 
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Symbols denoting the quantities are as follows: 

CO = • frequency at which the waves in waveguide propagate; 
= electronic volume charge density; 



jS = propagation constant appearing in the expression exp j(cot-jSz); 
c = velocity of light; 

e/m = 1.76 x 10^^ coulomb/Kg, electron charge to mass ratio; 

Hq = axial static magnetic field strength; 

“6 

Pq = 1.25 X 10 'henry/m; 

Eq = 8.86 X 10“^^ farad/m; 

4 



a = radius of the cylindrical waveguide. 

The explicit solution of the biquadratic equation is: 



2 2 

r To /,.2 ,.2 ,.2 Vc n /'„2 ,. 2 ^ 2 , 2 ^, 2 






2(k^-k^-k^ ) 
P cy' 






py p 



2 2 
k k 
. cy p 

- 



, f „2 , 2.2 ^^2 k „2 , 2,1 

-k ) +4^ — ^ (k -k )j 



(5) 



cy 



The plus sign in front of the square root will be consistently used for the 

2 2 
calculation of u^^ and the negative sign for the calculation of u^. 

2 

It is evident that we may multiply the Equation (5) by a throughout 
2 

and then obtain u expressed in terms of ka, k a, k a, fia which is 

P cy ^ 

useful for representation on diagrams since these values are then pure 



5 



numbers . 



The following relationship holds between the function u and the 

2 10 
function x used by Suhl and Walker : 



X 



2 

1,2 



2 

(ka)^-(kpa)^ 



( 6 ) 



It appears that the reason for the use of this expression by Suhl and 

Walker was the desire for unified treatment of both ferrite and plasma 

filled waveguides. Apparently propagation was not expected fOr values 

of fred[uency lower than the value of plasma frequency. A numerical 

solution of the equations of Suhl and Walker would have discovered the 

2 

plasmaguide modes that were later found by Trivelpiece and Gould . 

14 

Camus and Le Mezec use the notation; 




(7) 



and 



Z = 




(8) 



When boundary conditions are introduced it is best to normalize all 
frequencies to the frequency defined by; 



to 



0 



= c/a 



(9) 
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This means that quantitites k a , k , ka can be considered essentially 

P cy 

as to , to /to normalized with respect to to„. Whenever boundary 
p cy 0 

conditions are introduced the radius of the waveguide appears in the 



expression and both plasma frequency tOp and cyclotron frequency to^^ 

are specified by k a and k a respectively, i.e., in terms of the 

P cy 

normalizing frequency tOg = c/a. The solution of the boundary condition 

problem depends on the absolute value of to and to . However, when 

P cy 

2 2 

we consider only the quantities u^ and U 2 as functions of the frequency 
to and wave number ^ , i.e. , as functions in the tO“/S plane it is 
possible to normalize all frequencies to the plasma frequency (or alterna- 



tlvely to the cyclotron frequency) and the ratio of the cyclotron frequency 
to the plasma frequency becomes the only parameter that is left to be 
varied. This ratio is denoted in the following way: 



M = 



to 



cy 



co_ 



V’ 



( 10 ) 



Likewise, it is convenient to normalize the other quantitites with respect 
to the plasma frequency: 



Y = 



(ka) 

(kpa) 




( 11 ) 



X = 



(kpa) 




( 12 ) 
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2 2 

In this normalization the Hahn's equation for and U 2 reads: 



u 



ix.2 



(ka)' 

P 



— i _/.[2 (y2_i_m 2)+ M ](x^-Y>2CVr^-l) 






2 4X^Y^ 



(Y^-1) } 



(13) 



2 2 

Note that and U 2 are multivalued at the origin (X=0, Y=0). Its 

value depends on the path in the Y-X plane along which we approach the 

2 

origin*. It is seen without difficulty that the values of U are real 

X t L 

for values of X and Y exterior to the region in which the expression 
under the square root becomes negative, viz, , 



0 < Y < 1 



(14) 



^[/l+M^-Y^ - / 1-Y^ ]<X< -[Vi+M^-Y^ + \/ 1-Y^ J (15) 



2 2 

Within this region the values of and U 2 are complex conjugates. 
The boundary of the complex region is a closed curve touching the line 



Y = 1 at X = 1 and the line X = 



1+M^ ^ ^ /I + M‘ - ,, , 

at Y = / — .As M is 



M" V 2+M 

increased the complex region becomes narrower. This boundary is also 



shown on Figures 1 to 4 and 6 to 15. On the above boundary curve 
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f » 






f 



• W " 



I 






m* • 



i‘li 




i I 



I 







. I 






2 

Fig= 1 Contour curves for constant values of Re(U ) in the complex 
region of the co = plane at the value of M = 0o8o 
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2 

Fig, 2 Contour curves for constant values of Im(U ) in the complex 
region of the u.'' ~ S plane at the value of M = 0,8, 



10 



# 




2 

Fig. 3 Contour curves for constant values of Re(U ) in the complex 
region of the oo - $ plane at the value of M = 1, 6. 
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Fig. 4 Contour curves for constant values of Im(U ) in the complex 
region of the o) - j8 plane at the value of M = 1 . 6. 
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(Eqs. 14 and 15) the expression under the square root vanishes and = 
2 2 2 

U„. The values of U. = U_ on the boundary curve are: 



on 




[^1 + M^-Y^ 




(16) 



M 



2 , 2 2 ^ 
2Y (1+M -Y ) 



-l] /(M^+1-Y^)(l-Y^) } (17) 



and on 






y 



2 2 
1+M -Y 






Z 1 



(18) 






2Y^(l+M^-Y^) 



-1,] \/(M^+1-Y^)(1-Y^) } (19) 



If the origin (X=0, Y=0) is approached along the curve defined in Eq. 16 
2 

the value of U at zero is given by 

; ( 20 ) 

y 1+M^ 

if on the other hand the origin is approached along the curve defined in 

2 

Eq. 18 the value of U at zero is given by 





\ 

1+M^ 



( 21 ) 
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The equations of contour lines Re(u ) = const and Im(ir) = const within 

2 

the region defined by Eq. 14 and 15 are readily found. For Re(U ) = constant 



X = 




[ -2(M^+1-Y^) ] 



( 22 ) 



L ^2 



for Im(U ) = constant 




(M^+l-Y^)+(l-Y^)+2 




(l-Y^)(l+M^-Y^)-(M^+l-Y^)^(Im(U^) (23) 



The curves given by Eq. 22 are single valued functions of Y while the 

curves given by Eq. 23 are double valued; they are closed curves 

similar to the boundary of the region. Figures 1 to 4 show representative 

contour diagrams inside the complex region. It is further of interest to 

determine whether and in which points on the boundary curve the value 
2 2 

of U vanishes. Evidently U =0 at Y=1 and X=l, other points are 

found by setting the expressions of Eq. 17 or 19 equal to zero, and 

2 

arrange the terms in powers of Y . Thus we obtain; 



The roots of this equation are real as would become apparent if Eq. 24 were 
solved; it is possible to obtain them by solving them numerically. However, 




(24) 



14 




^ .-Sr * •« : 



t 






•« • # ' • 



« 












a more elegant way to study its roots is available. It will be recalled that 
2 2 • 

either or U 2 can vanish on the following two curves; 




These two curves lie outside of the complex region. Equation 25 re- 
presents a curve with a discontinuity at Y = M, its upper branch inter- 



cepts the Y axis (corresponding to co axis) at Y= + ~ 




+ 1 ; 



Eq. 26 represents a continuous curve with Y 



Y = - ^ + y*— + 1. Values of and U 

Equations 25 and 26 can be readily evaluated. 
Eq. 25 we have: 



intercept at 
> * 

I on curves given by 
On the curve given by 



y(Y^-1-M^) 

Care must be exercised in taking the square root in this expression. 

We shall here take the square root as the opposite operation of squaring. 
The expression under the square root is always positive but the expression 
squared can be either positive or negative. For the sign of the square 
root we take the sign of the squared quantity which we must first evaluate. 



| - Y^-MY^-Y+M/^ ^ Jq Y^-MY^-Y+M^ yj 



3 ._.2 

Y-M 



(27) 



15 




Fig. 5 Plot of the parameter M as a function of Y(=a’/top) illustrating 
the solutions of the cubic of Eg, 28. 
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The roots of the cubic equation 



Y^-MY^-Y+ ^ = 0 



may be determined as follows. From Eq. 28 we have 



(28) 



M = 



Y(Y-1)(Y+1) 

sfi 'JT 



(29) 



The value of M can be readily plotted as a function of Y as is shown 
on p'ig. 5 . From Fig. 5 it is seen that one root of Eq. 28 de- 

signated as first root, is always negative and that the third root is 
always found between the following limits: 

. If we limit our interest to the first quadrant of the o)-/3 plane 
(Y >0, X > 0) and to real values of X these two roots need not be 
considered. The second root, on the other hand, lies between the values 

0 < Y, < i (31) 

This means that a point exists on the lower branch of the curve given by 

Eq. 25 with an ordinate Y which satisfies the cubic of Eq. 28 and as a 

2 2 

consequence both and are equal to zero in that point. If both 
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and U 2 are equal to zero then they are also equal to each other and 

consequently lie on the curve given by £q. 18 « i.e. < on the right 

2 2 

boundary of the region of complex values of and U 2 . Hence the 



curve of Eq. 25 has one point in common with the right boundary of the 
complex region given by Eq. 18 and the two curves are tangent to each 
other in that point. We shall call this point the lower tangent point. Now 
for values of Y between zero and the second root Y 2 of Eq. 28 the 
squared expression in Eq. 27 is negative, viz. , 



‘32 M 
Y -MY -Y+ -~ > 0 



(32) 



Y - M < 0 



(33) 



if 



0 < Y < Y, 



2 



(34) 




12 2 
Y(Y -I-m"^) 




Y^-MY-YHNy2 

Y-M 



)}=0 



(35) 



and for U^: 



2M(Y^-MY-Y+ M/ 2) 
^ Y(Y^-1-M^)(Y-M) 



(36) 
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41 




When, on the other hand, the value of Y is greater than the value of the 



laoond root tho following relations hold: 

Y^-MY^-Y+ ^ .< 0 (37) 

Y-M < 0 (38) 



Y^< Y < M 

The squared quantity in Eq. 27 is then positive and we have 

^ 2M(Y^-MY^-Y+M/2) 

^ Y(Y^-1-M^)(Y-M) 






M 



Y(Y^-1-M^) 



{ 



Y^-MY~Y+M/2 Y^-MY^-Y+M/2 'i 

Y-M " Y-M J" ^ 



(39) 



(40) 



(41) 



Thus it is seen that on the curve given by Eq. 35 between the origin 

(X=0 , Y=0) and the point in which this curve touches the boundary of the 

complex region the function is zero and the function has a value 

given by Eq. 40, on the rest of the curve, beyond the lower tangent point, 

2 

the situation is reversed and the function U 2 is zero while the function 
2 

is given by Eq. 40. Likewise it can be seen that on the upper branch 

2 

of the curve given by Eq. 25 the value of the function U 2 is always zero 

2 

and the value of the function U 2 is given by Eq. 40 for positive values 
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of y. 

2 2 

In the same manner we obtain the values of and U 2 on the 
curve given by Eq. 26. 



2 M r Y^+MY^-Y-M^ A Y^+MY^-Y-M^ ^2 

1,2 y(Y^-1-M^) ^ ' ^ 



( 42 ) 



To determine the sign of the square root it is again necessary to study the 
behavior of the cubic polynomial and the roots of the equation 



y3*My2-Y- 0 



(43) 



It is immediately obvious that this equation becomes identical to Eq. 28 
if the variable Y is replaced by -Y. The roots are then as shown on 
Fig. S only with the positive and negative Y axis interchanged. 

From the same considerations as before only the positive root, designated 
henceforth as the third root (of Eq. 43) is of interest. The third root- 
lies between the values 

^<Y3<1 '(44) 

Thi^ again, means that a point exists on the curve given by Eq. 26 with 

an ordiante Y^ which satisfies the cubic of Eq. 43 and consequently both 
2 2 

and are equal to zero in that point and, according to the argument 
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-rT’' 



presented before, this point, too, lies on the curve given by Eq. 16, i.e. , 



on the left boundary of the region of complex values of and U 2 < Thus 
the curve of Eq. 26 has a point in common with the left boundary of the 
complex region given by Eq. 16, the two curves are tangent to each 
other at that point. We shall call this point the upper tangent point. Now 
for values of Y between the ordinate intercept of the curve given by the 
Eq. 26 and the third root, viz. , 






(45) 



the squared expression in Eq. 42 is negative 



3 2m 
Y +MY -Y- -f- < 0 



(46) 



From where follows for U, 



-2M r Y^+MY^-Y-M/2 -I 



(47) 



a nd for U„ 






M 



Y(Y^-1-M^) 



[■ 



_ Y^+MY^-Y-M/ 2 Y^+MY^-Y-M/2 > n 

Y+M " v' Y+M J J 



(48) 
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when, on the other hand, the value of Y is greater than the value of the 



third root the following relations hold 

+ MY^ - Y - > 0 (49) 

Y > Yg . (50) 



The squared quantity in Eq. 42 is then positive and we have for 



U?= 



^ Y(Y^-1-M^ 



L” Y+M V Y+M J J 



= 0 



(51) 



and for U,. 



2^_^2M r 

^2 2 2 
Y(Y -1-M ) 



Y^+MY^-Y-M^ 1 



Y+M 






(52) 



In the same manner it can be shown that on the left boundary of the 

2 

complex region given by Eq. 16 the values of U are negative for values 

of Y samller than the value of the third root (Y< Y^) and positive for 

values of Y greater than the value of the third root (Y^ < Y < 1) and 

2 

that on the right boundary given by Eq. 18 the values of U are positive 
for values of Y smaller than the value of the second root (Y < and 
positive the values of Y greater than the value of the second root 
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(Y^ < Y < 1). Moreover it can be seen that on the line Y=1 the values 
2 2 

of functions and U 2 are as follows: 

In the interval: 

Y=1 , 0 < X < 1 (53) 




2 2 
U2 = l-X 



and in the interval: 

Y=l, KX < ® 



(54) 

(55) 

(56) 



.= 1 - X^ (57) 

= 0 (58) 



To show that the roots of Equations 28 and 43 are also the roots of 

2 2 

Eq. 24 obtained for the condition that both and U 2 vanish we 

multiply the cubic polynomials of Equations 28 and 43; the resulting 

/ 

bicubic polynomial is Just that of Eq. 24. 

2 2 

Since in the complex region the functions and U 2 are 

complex conjugates they can be represented by their modulus p and 
argument <p , viz . ; ^i~p / ^2~P / * Note that we have 

here and U 2 instead of and U^. Contour curves of constant 
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modulus and argument of the functions and U 2 are of greater 

practical value than the contour curves for constant real and imaginary 

2 16 
part of U . Tables of Bessel functions of complex arguments give the 

argument of the Bessel functions in polar form. 

jr' 

The modulus p is readily found from Eq. 13; 




1 

, 2 2.2 

4(Y -1-M ) 




M 



)+ ^)(x''^-yV2(Y^-1)]' 



.2 „2 



(59) 



I 




(l-Y^)-(X^-Y^)^ ] } 



To find the equation of a contour curve in the Y-X (a? - jS ) plane for 

constant value of p the above equation is expanded in powers of X 

2 

and a quadratic equation in X is obtained. 



c/ + c/ + c„ - p'* = 0 



(60) 



The coefficients C^, depend only on Y and M and are given 
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by the expressions; 











( 61 ) 





2 2 2 

Aix -1 -mY 






(63) 



From Eq. 60 we obtain the explicit form of’ the equation for the contour 
curve of constant modulus , p . 




Obviously X must be real and positive. Some representative contour 
curves given by the above equation are shown on Figures 6 and 8. 
The argument <p is also found from Eq. 13: 



M 

,.2 



tan (2(p )= 



4 

M 



-[2(Y^-1-M^)+ ]- 2(Y^-1) 



(65) 



To find the equation of a contour curve for constant argument in the Y'-X 
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(o)“j3) plane we use the notation T = tan {2ip) and expand the Eq. 65 in 

2 2 
powers of X so that a quadratic equation in X is obtained. 






( 66 ) 



The coefficients Cq, C^/ depend only on Y, M, and T and are 
given by the following expressions: 



Cq=tV[ 2(Y^-1-M^)+ ^]^-4T^(Y^-l)Y^j2(Y^-l-M^)+ ^ ]+4T^(Y^-1 )^+m'^ (67) 



C2- -2Y^T^[2(Y^-1-M^)+ ^]^+4T^(Y^-l)[2(Y^-l-M^)+ *^^(1-Y^) (68) 



2 — tv/ 

M 



C^= T^[2(Y^-l-M^) + ^ ]^ + 



(69) 



From Eq. 66 we obtain the explicit form of the equation for the contour 
curve for constant argument. 




Again, X must be real and positive. Some representative contour curves 
given by the above equation are shown on Figures 7 and 9. 
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Fig. 6 Contour curves for constant values of the modulus o of the function 
U = o ^ in the complex region of the a' - plane at the value of 
M=0.8. 
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Fig. 7 Contour curves for constant values of the argument o of the 
function U = o ^ in the complex region of the a' - ^ plane 
at the value of M = 0.8, 
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M 

(t)o 



Fig. 8 Contour curves for constant values of the modulus o of the 
function U = ^ in the complex region of the Lc - ^ plane 
at the value of M = 1 . 6 . 
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Fig. 9 Contour curves for constant values of the argument o of the 
function U = oZl in the complex region of the a? - jS 
plane at the value of M = 1 , 6. 
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It is interesting to note the curves delineating various zones in 



the complex region. On tie left boundary of the complex region between 
th« origin and the upper tangent point the value of the argument of U !■ 
7 t /2, in the upper tangent point the value of the argument is not defined 
and can assume any value, and between the upper tangent point and the 
point (X=l, Y=l) the value of the argument is zero. In this point, again, 
the argument is not defined. 

On the right boundary of the complex region between the origin 

and the lower tangent point the value of the argument is zero, in the 

lower ‘tangent point its value is not defined and between the lower tangent 

point and the point (X=l, Y=l) the value of the argument is ir/2. 

2 

On curves where Re(U ) = 0 the value of the argument is 7t/4. 

Contour curves converge in the upper and lower tangent points, in the 

point (X=l , Y=l), and in the origin. 

It is now possible to obtain a picture of the behavior of the 
2 2 

functions and U 2 in various regions of the Y-X (to - jS) plane. 

The curves given by Equations 14, 15. 25, 26 and the lines Y=1 and' 



whether the value of M is less or greater than unity. It has already 

been pointed out that within the region defined by Equations 14 and 15 
2 2 

the functions and U 2 assume complex conjugate values . In the 

2 

other regions the following has been found concerning the sign of 




divide the Y-X (o) - j3) plane into 12 or 13 zones depending 



and U^: 
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Fig» 10 Zones in the w - ]8 plane in which the functions and ul 

are of the same or of opposite sign. General behavior for 
values of M < lo 
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Fig. 11 Zones in the w- $ plane in which the functions 
2 2 

Uj and U2 are of the same or of opposite sign. 
General behavior for values of M > 1. 
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At Values of M smaller than unity, in the region above the upper branch 

2 2 

of the curve given by Eq. 25. Both and U2 are positive. Between 
the upper branch of the eurve given by £qi 25 t and the ourve given by 

/ 2 2 2 

Eq. 26, above the straight line Y= l+M ; > 0, U2 <0. On the 

/ 2 2 2 
line Y= v 1+M the function U2 has a pole while the function 

is continuous. To the right of the curve of Eq. 26 and above the line 

j 2 2 2 

Y=^l+M both and U2 are negative. Between the lines 
Y=1 and Y= J 1+M^ both and are positive to the left of the 

2 

curve of Eq. 26 while to the right of the curve of Eq. 26 the value of 
* 2 

is negative and the value of U2 positive. Under the line Y=1 and 

2 

above the curve given by Eq. 26 the value of is negative while that 
2 

of U2 is positive. Under the line Y=1 and under the curve of Eq. 26 

and above the left boundary of the complex region given by Eq. 16 both 
2 2 

Uj^ and U2 are positive. Under the curve given by Eq. 26 and to the 

2 

left of the boundary of this complex region given by Eq. 16 both U^^ and 
2 

U2 are negative. Between the right boundary of the complex region 

given by Eq. 18 and the lower branch of the curve of Eq. 25 both U^ and 
2 

U2 are positive. Under the lower branch of the curve of Eq. 25 the function 
2 2 

Uj^ is negative while U2 is positive. To the right of the right boundary 

of the complex region given by Eq.l8 under the line Y=1 and above the 

2 2 

lower branch of the curve of Eq. 25 both U^ and U2 are negative. 

Figure 10 shows the zones and signs of the functions U^ and U^ for 
the case when M < 1. 
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For values of M greater than unity the lower branch of the curve given 
by Eq. 25 and the line Y=1 intersect and modify the picutre with the 

2 

addition of another zone. In such cases the behavior of functions Uj^ 

2 

and is modified as follows. To the right of the lower branch of the 

» 

2 

curve of Eq. 25 and below the line Y=1 the function is negative 
2 

while U 2 is positive. In the newly formed triangular zone bounded by 

the right boundary of the complex region given by Eq. 18, the line Y=l, 

2 2 

and the lower branch of the curve given by Eq. 25 both and U 2 are 
negative, likewise they are negative in the region above the line Y=1 

and bplow the lower branch of the curve given by Eq. 25. In the region 

✓ 

to the right of the curve of Eq. 26, above the line Y=1 and above the 



lower branch of the curve of Eq. 25 and below the line Y=v M +1 

2 2 
the function is negative while the function. U 2 is positive. The 

2 2 

behavior of and in the other zones remains the same as for 

cases when M < 1. Figure 11 shows the zones and signs of functions 

2 2 

and U 2 for this case (M> 1). 

The partition of the to -/S (Y-X) plane into zones as described above 

is quite general and is valid for all possible cases. These cases fall • 
under one of the two types, viz. , those with M< 1 and those with 
M> 1. There is, however, some variation in the sequence of special 
frequencies Y^, ~ ^ ■*" \/ \/ / etc.) depending in 

which interval lies the value of M. The sequential order of these 
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special frequencies depend only on the values of M. The possible cases 
are listed below. 



0< M<^ 

0 < < M < 




< 





(71) 

(72) 



,^ < M < 1 (73) 

V2 







<Y^ < M < 1 





(74) 



1 < M < 1.5 . (75) 

° ^^2^ “ T +1 < Y3 < l.< M <yM^+l < ^ + +1 (76) 



1.5 < M 


/ 








(77) 






L < M 


/m^+i <^ + ^ 


14 


(78) 






In addition to the general behavior of the functions U., and U„ or, 

♦i i ^ 

if roots are taken, and U^, it is also desirable to have contour 
representations of these functions in the to ~ j8 (Y-X) plane where 
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and have either real or imaginary values, i.e. , outside of the complex 



region. 



The equation of contour curves in the co ~ /S (Y-X) plane for which 
U is a real or an imaginary constant can be readily obtained by solving 
Eq . 3 for X . 



This is a simple solution of a biquadratic equation in X. For contour 



curves in the first quadrant of the ui - ^ (Y-X) plane only positive 
2 

values of X are acceptable. There appears to be no simple rule to 
determine whether the curves obtained from Eq. 79 are the contour 
curves of or U 2 . This question can be resolved by calculating with 
the obtained values of Y and X the desired function or \J^, as the 
case may be, using the rule concerning the selection of the sign as 
stated above and then determining whether the calculated value is equal 
to the value of U originally set in Eq. 79.. If the calculated value of 
Uj or U 2 is equal to the value of U set in Eq.79 the point in question 
lies on the curve = constant or U 2 = constant respectively. 

Equation 3 can also be solved with respect to Y albeit not so readily 
since the arrangement in powers of Y leads to a biquartic eqxiation. 

This biquartic equation can have complex roots and the selection of proper 
curves is even more laborious in this case. There appears to be no point 




(79) 
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in solving the Eq. 3 for Y, therefore, since the solution with respect to 
X can be obtained by simple algebra. 

Some representative contour curves for various real and imaginary 
values of and U 2 are shown on Figures 12, 13, 14, and 15. 



4 
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Fig. 12 Contour curves for constant real and imaginary values of U, 
in the co - /3 plane at the value of M = 0.8. 



I 
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Fig. 13 Contour curves for constant real ana imaginary values oi U, 
in the “ B plane at the value of M = 0.8. 
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Fig. 14 Contour curves for constant real and imaginary values of 
in the a' - 8 plane at the value of M = 1.6. 
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Fig. 15 Contour curves for constant real and imaginary values of U 2 
in the co - S plane at the value of M = 1.6, 

41 
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